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Abstract 

The approach to the constructing explicit solutions of the recurrence relations for 
multi-loop integrals are suggested. The resulting formulas demonstrate a high effi- 
ciency, at least for 3-Toop vacuum integrals case. They also produce a new type of 
recurrence relations over the space-time dimension. 



1 Vacuum case 



Recently [1,2] a new approach to implement recurrence relations [3] for the 
Feynman integrals was proposed. In this work we extend the general formulas 
for the solutions of the recurrence relations to the multi-loop case. Let us 
consider first vacuum L-loop integrals with N = L(L+ l)/2 denominators (so 
that one can express through them any scalar product of loop momenta) of 
arbitrary degrees: 
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The recurrence relations that result from integration by parts, by letting 
(d/dpi) ■ pk act on the integrand [3], are: 
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where lfB(. . . , n c , . . .) = B{. . . , n c ± 1, . . .), in particular l^n a = n a ± 6°. 
Using the relations 
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The common way of using these relations is step-by-step reexpression of the 
integral (1) with some values of rii through a set of integrals with shifted values 
of rii, with the final goal to reduce this set to a linear combination of several 
"master" integrals N k (D) with some "coefficient functions" F k (n,D): 



B(n, L>) = £ F k (n, D)N k (D) 



Nevertheless, to find proper combinations of these relations and a proper se- 
quence of its use is the matter of art even for the tree-loop integrals with one 
mass [4]. Then, even in cases when such procedures were constructed, they lead 
to very time and memory consuming calculation because of large reproduction 
rate at every recursion step. Instead, let us construct the F k (n,D) directly 
as solutions of the given recurrence relations. Note, that if we find any set of 
the solutions, we could construct F k (n,D) as their linear combinations. Let 
us try the solution of (3) in the following form: 
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where integral symbols denote iV subsequent complex integrations with con- 
tours which will be described later. Acting by some operator n a I a+ ) (all 
decreasing operators should be placed to the left) on (4) and performing the 
integration by parts one gets (s.t. are surface terms): 
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^Z7T I J J J X-y ' 'Xjy 

So, if we choose the g(x a ) as the solution of Oi(x a ,d a )g(x a ) = and cancel 
the surface terms by proper choosing of integration contours (for example, 
closed or ended in the zero points) we find that (4) is a solution of relations 
Oi(I~, n a I a+ )f k (n) = 0, and different choices of contours correspond to differ- 
ent solutions. 

The differential equations for (3) have the solution g(x a ) = P(x a +fi a ) lyD ~ L ^ 1 ^ 2 , 
where 

N 

P(x a ) = det(J2(A-% l)Xa ) 

a=l 

is the polynomial in x a of degree L, so we get the desirable solutions of (3): 

Finally, let us derive from (3) the recurrence relations with D-shifts. Note 
that if f k (iii, D) is a solution of (3), then by direct substitution to (3) one can 
check that P(I~ + fia)f k {ni, D — 2) also is a solution. Hence, if f k {n i) D) is a 
complete set of solutions, then 

f k (n u D)=J2 S k (D)P(r + ^)f n (n t , D — 2), 

n 

where the coefficients of mixing matrix S is numbers, that is do not act on 
rij. For the solutions (5) the S is the unit matrix (the increasing of D by 2 
leads to appearing of factor P(x a ) in the integrand of (5)), but the desire to 
come to some specific set of master integrals may lead to nontrivial mixing. 
These relations look different from recently proposed in [5] , although further 
investigations can give some connections with them. 

To check the efficiency of this approach we evaluated (using REDUCE) the 
first 5 moments in the small q 2 expansion of the 3-loop QED photon vacuum 
polarization. The 3-loop contribution to the moments are expressed through 
about 10 5 three-loop scalar vacuum integrals with four massive and two mass- 
less lines. The integral (5) in this case can be solved to finite sums of the 
Pochhammer's symbols (see [1]). Moreover, it is not necessary to evaluate 
these integrals separately. Instead, we evaluated a few integrals of (5) type, 



3 



but with P D / 2 2 producted by a long polynomial in Xi (the results see in [1,6], 
they are in agreement with QCD calculations [7] made by FORM). 

The comparison with the recursive approach shows a reasonable progress: the 
common way used in [6] demands several CPU hours on DEC-Alpha to calcu- 
late full D dependence of the first moment, and further calculations became 
possible only after truncation in (D/2 — 2). In the present approach the full 
D calculation for the first moment demands a few minutes on PC. 



2 Non-vacuum case 



Suppose that integrals (1) depend on R external momenta pi (L < i < L + R). 
The number of the denomenators are now Ni = L(L + l)/2 + LR, and the 
number of additional ("external") invariants are N 2 = R(R + l)/2. Let us 
expand the integrals in formal seria over "denominator-like" objects D a of (2) 
type with a — Ni + 1, .., N± + N 2 , depending on external momenta only: 
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We define such general expansion in order to write the recurrence relations 
in compact form, in practice the coefficients and fx a may be very simple. 
The expansion with negative rii corresponds to the large momenta expansion, 
with positive ones to the expansion near points /i a m 2 . The rii can also be 
noninteger, but with unit shifts. 

Acting by (d/dpi) ■ pk, (i = 1, . . . , L; k = 1, . . . , L + R) on the integrand 
we get N\ recurrence relations. The additional N2 relations we get acting by 
Pk ■ (d/dpi), (i, k — L + 1, . . . , L + R) on both sides of (6). These new relations 
look like the old ones with only exception that they have no terms proportional 
to space-time dimention D. The complete set of recurrence relations is now 

((D-L-R-l)S i k -(D-R-l)S i k )b(n,D) = 

N!+N 2 L+R 
a,d=l 1=1 
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where 5 k ={5 k if i, k > L, else 0). The corresponding differential equations have 
the solution g(x a ) = g'(x a + fi a ), where 

g'(x) = det((A-% l) x a )'^^det ((A-% l) x a y £ ^ 1 } (7) 

and det denotes the minor with k, I > L. So, one can use the representa- 
tion (4), but the problem of resolving it to explicit formulas demands futher 
investigations. 

Finally note that one can formally obtain the formulas (5, 7) by "change of 
integration variables" from loop momenta to " denomenator-like objects" D a . 
The weight function for this change is 

J d D Pl ■ ■d D p L Y[5(D l /m 2 - Xi) (xdet{{A- 1 ) a (kl) {x a + fia)) 2 ^ 1 . 
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